m ■ m m 

!MA>»A 

Technical 

Paper 

2758 


September 1 987 


a i „r t- i i 

/^ppiI^dUUIl <J1 lUlUUiCIltC 

Modeling to Predict 
Surface Heat Transfer in 
Sta.pnation Flow Region 

o o 


W 1 


f ^ ir^nl 


i inr 

VJll C/UiUi 


r\,i 

y rrx lvjls^x 

j 


Chi R. Wang and 
Frederick C. Yeh 


1% I 

iu/o/\ 


NASA 

Technical 

Paper 

2758 

1987 


Application of Turbulence 
Modeling to Predict 
Surface Heat Transfer in 
Stagnation Flow Region 
of Circular Cylinder 


Chi R. Wang and 
Frederick C. Yeh 

Lewis Research Center 
Cleveland, Ohio 


NASA 

National Aeronautics 
and Space Administration 


Scientific and Technical 
Information Office 


Summary 

A theoretical analysis and numerical calculations for the 
turbulent flow field and for the effect of free-stream turbulence 
on the surface heat transfer rate of a stagnation flow are 
presented. The flow field considered is the region near the 
forward stagnation point of a circular cylinder in a uniform 
turbulent mean flow. The emphasis of the present study is on 
the modeling of turbulence and its augmentation of the surface 
heat transfer rate. The free stream is steady and incompressible 
with a Reynolds number of the order of 10 5 and a turbulence 
intensity of less than 5 percent. In the analysis, the flow field 
is divided into three regions: (1) a region significantly away 
from the cylinder where the mean velocity is uniform and the 
turbulence is homogeneous and isotropic, (2) an external 
inviscid flow region with a constant rate of the mean velocity 
variation where the turbulence is distorted by the mean flow 
velocity, and (3) an anisotropic turbulent boundary layer region 
over the cylinder surface. The turbulence modeling techniques 
used are the k-t two-equation model in the external flow region 
and the time-averaged turbulence transport equation in the 
boundary layer region. The turbulence double correlations, 
the mean velocity, and the mean temperature are solved 
numerically from sets of finite difference equations. These 
finite difference equations are derived from the flow-field 
governing equations. The solution is obtained from matching 
the turbulence kinetic energy along the boundaries between 
the flow regions. The results and their comparisons with the 
experiments show that (1) depending on the free-stream 
turbulence characteristics, the turbulence kinetic energy is 
amplified or attenuated along the stagnation-point streamline, 

(2) the free-stream turbulence penetrates into the flow region 
near the surface and induces a high surface heat transfer rate, 

(3) the turbulence kinetic energy increases continuously in the 
mean flow direction along the boundary layer edge and, near 
the stagnation point, the rate of variation in the Reynolds 
normal stress is different for each component, (4) depending 
on the level of the turbulence kinetic energy, different surface 
heat transfer rates may occur near the stagnation point, and 
(5) the analysis provides a procedure in computing the surface 
heat transfer rate near a stagnation point as functions of the 
turbulence longitudinal microlength scale, the turbulence 
intensity, and the free-stream Reynolds number. 


Introduction 

The momentum and thermal flow fields near the forward 
stagnation point of a circular cylinder in turbulent flow have 
been the focus of considerable study (refs. 1 to 12). One 
objective has been to understand the enhancement of the heat 
transfer due to the free-stream turbulence near the leading edge 
of an airfoil. Existing experiments (refs. 3 and 4) have shown 
that when the turbulence intensity is 3 to 5 percent the local. 
Nusselt number increases by 75 percent as compared with 
laminar flow conditions. In the initial stage of turbine blade 
design, the blade leading edge is often modeled as a circular 
cylinder in a turbulent free stream. With the current emphasis 
on high turbine inlet temperature, it has become even more 
important to understand the effect of turbulence on the surface 
heat transfer, particularly at the leading-edge stagnation region. 
In fact, prediction of the stagnation-region surface heat transfer 
becomes vital to the turbine blade cooling design. To estimate 
this surface heat transfer rate, correlations derived from 
experimental data may be used. Detailed experimental data 
can be found in the paper by Lowery and Vachon (ref. 5). 
The alternative is to rely on the analysis to obtain the surface 
heat transfer rate. The objective of the present study is to 
develop an analytical method that can predict the effect of the 
turbulence on the surface heat transfer rate. 

Smith and Kuethe (ref. 6) used the two-dimensional 
boundary layer theory to analyze the flow at the stagnation 
point on a circular cylinder. The eddy viscosity was assumed 
to be proportional to the turbulence in the free stream. 
Reynolds’ analogy was used. This analysis indicates a large 
effect of turbulence on heat transfer. A similar analytical 
approach was used by Gorla (ref. 7) to investigate the effects 
of the unsteadiness in the free stream and its turbulence on 
heat transfer at a stagnation point. Other studies (refs. 8 to 
11) show that turbulence amplification, because of vortex 
stretching, can also influence the thin boundary layer and, 
therefore, the surface heat transfer rate. 

To model the turbulence effect, Hijikata et al. (ref. 12) 
performed a theoretical and experimental study of the 
stagnation-point anisotropic turbulence. An equation for the 
anisotropic turbulence was added to the k-e two-equation 
turbulence modeling to analyze the measurements. Their 
analysis shows that the turbulence length scale affects the 
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surface heat transfer rate. Traci and Wilcox (ref. 13) 
analytically studied the turbulence effect on the stagnation- 
point flow field using Saffman turbulence modeling (ref. 14). 
The mean velocity gradients within the inviscid flow were 
found to be important in analyzing the stagnation-point heat 
transfer. A difficulty in this type of analysis is attaining a 
proper match of the boundary conditions between the different 
flow regions. Strahle (ref. 15) tried to resolve analytically the 
turbulence matching conditions. Analytical relations between 
the turbulence kinetic energy and its dissipation rate at the 
stagnation point and in the free stream were found. On the 
basis of the method of invariant modeling (refs. 16 and 17), 
Wang (ref. 18) proposed a theoretical boundary layer flow 
analysis and numerical computational procedure to investigate 
the turbulence, momentum, and thermal fields within the 
stagnation flow region. Ad hoc turbulence boundary conditions 
were imposed along the boundary layer edge. The numerical 
calculation showed that rapid amplification of the Reynolds 
normal stresses along the boundary layer edge was required 
to predict the approximate surface heat transfer rate as reported 
in existing experiments. 

From the experimental results and analytical methods of 
these existing studies, it is possible to pursue a detailed 
theoretical and numerical analysis of the forward stagnation 
flow field of a circular cylinder in a turbulent free stream. 
Reported herein is such an analysis, which provides a 
procedure to compute the surface heat transfer rate as functions 
of the free-stream turbulence longitudinal microlength scale, 
the turbulence intensity, and the Reynolds number. The 
emphasis is on the modeling of turbulence and its augmentation 
of the surface heat transfer rate. 

The stagnation flow field is divided into three distinct 
regions: (1) a region away from the cylinder where the mean 
velocity is uniform and the turbulence is assumed to be 
homogeneous and isotropic, (2) an external inviscid flow 
region which has a constant rate of mean velocity variation 
and the mean velocity distorts the turbulence, and (3) an 
anisotropic turbulent boundary layer region over the cylinder 
surface. For turbulence modeling, the k-e two-equation 
modeling and the time-averaged turbulence transport equation 
are used respectively in flow regions 2 and 3. These turbulence 
modeling equations and the mean flow conservation equations 
are used to formulate the theoretical analysis. The turbulence 
double correlations, the mean velocity, and the mean 
temperature are solved numerically from sets of finite 
difference equations. The solution is obtained from matching 
the turbulence kinetic energy along the boundary between the 
flow regions. These analytical results are compared with the 
measurements from existing experiments. 

Symbols 

A a coefficients of equation (Cl), i = 1 to 9 
a mean velocity parameter 


£, coefficients of equation (C6), i = 1 to 4 
c'i empirical constants, i = 1 to 5 
c p specific heat 

Cp empirical constant 

D cylinder diameter 

E dimensionless turbulence dissipation rate 

F dimensionless mean flow streamwise velocity 

f parametrical functions of <p, i = 1 to 3 

G flow variable of equation (Cl) 

g turbulence longitudinal correlation 

H dimensionless mean temperature 

K dimensionless turbulence kinetic energy 

k turbulence kinetic energy 

( spatial distance 

Nu Nusselt number 

P mean static pressure 

Pr laminar Prandtl number, 0.7 

p pressure fluctuation 

R sum of dimensionless Reynolds normal stresses, 

R x + R y + R z 

R x dimensionless Reynolds normal stress along x direction 

R y dimensionless Reynolds normal stress along y direction 

R z dimensionless Reynolds normal stress along z direction 

Re Reynolds number 

S dimensionless Reynolds shear stress 

T mean temperature 

Tu free-stream turbulence intensity, 

t temperature fluctuation 

U,V mean velocity components 

V r transformed velocity of V (eq. (13)) 

v, turbulent dynamic viscosity 

m,v,w velocity fluctuations 
x,y,z physical coordinate system 

5 stagnation-point laminar boundary layer thickness 
e turbulence dissipation rate 

tj transformed coordinate along y direction 

9 dimensionless turbulence correlation, — vt/U e T e 
k thermal conductivity 

A turbulence modeling length scale 

X turbulence longitudinal microlength scale 

v kinematic viscosity 

£ transformed coordinate along x direction 

p density 

<p azimuth angle 
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Subscripts: 


v F , t f . Re 


e external or boundary layer edge condition 
F free-stream condition 

m grid point along £ or jj direction 

o,e boundary layer edge condition at initial station 

w surface condition 

0 condition near stagnation point 

1 condition at location 1 

2 condition at location 2 of boundary layer edge 

3 condition at location 3 of boundary layer edge 

Superscript: 

— time mean average 


Analytical Formulation of Flow Regions 

The flow region of interest is around the forward stagnation 
point of a circular cylinder in a turbulent free stream. Figure 1 
schematically represents the flow field. The diameter of the 
cylinder is D and its axis is perpendicular to the free-stream 
mean velocity V F . The free-stream Reynolds number is defined 
as Re = V F D/v, where v is the kinematic viscosity. The 
velocity fluctuation in the free stream is v F , and the free- 
stream turbulence intensity is defined as Tu = Vv J/V F . The 
free-stream turbulence longitudinal microlength scale is 
denoted by \ F . The cylinder is assumed to have a constant 
surface temperature T w . The free-stream temperature is T F . 

The stagnation point is taken as the coordinate-system origin; 
x and y denote the distances parallel to and normal to the 
cylinder surface; and z denotes the distance parallel to the 
cylinder axis. The mean flow field is assumed to be steady, 
incompressible, and two dimensional. The mean velocity 
components are U and V along the x and y coordinates. The 
turbulence is three dimensional. The velocity fluctuations are 
denoted by u 2 , v 2 , and w 2 along the x, y, and z directions. 
To formulate the theoretical analysis, the flow field is divided 
into three flow regions. The governing equations of the mean 
flow and the turbulence are derived for each flow region. 
Following the laminar flow analysis, the relations 

i * u 

— dx and t]{x,y) = — (1) 

o " ( 2£) 0 ' 5 

are used for coordinate transformation. The analytical 
formulation, together with the assumptions and the boundary 
conditions, for each flow region is described here. 



Figure 1.— Schematic of flow field of interest. 


Region 1 

This flow region is the uniform turbulent free stream away 
from the cylinder. The properties of the mean flow and the 
turbulence are not influenced by the existence of the cylinder. 
Appropriate conditions of the turbulence kinetic energy k and 
the turbulence dissipation rate e are required to characterize 
the turbulence. The turbulence is assumed to be homogeneous 
and isotropic. The free-stream turbulence kinetic energy and 
the turbulence dissipation rate are given by the following 
equations (ref. 19): 


k F — 



and 



( 2 ) 


The turbulence longitudinal microlength scale can be 
determined from the isotropic turbulence theory and two-point 
turbulence correlations. The two-point turbulence longitudinal 
correlation g, the spatial distance f, and the free-stream 
turbulence longitudinal microlength scale \ F are related by 
the following equation: 
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^)=‘-72 

A/r 


The boundary conditions are defined in the following: 

(1) At T7 = 77 ] , 


with small values of t 
Region 2 

In this region, the mean flow velocity is disturbed by the 
cylinder. The mean velocity variation distorts the turbulence. 
However, the turbulence does not influence the mean flow 
velocity. Thus, the mean flow is represented with the two- 
dimensional inviscid flow solution near the stagnation point. 
The mean flow velocity components are 

U = ax and V = - ay (4) 

with a = 4 V f /D. The turbulence kinetic energy variation 
along the stagnation-point streamline is of interest in this study. 
The turbulence is modeled with the k-e two-equation 
turbulence model. A summary of this modeling technique is 
given in appendix A. For convenience, the x direction diffusion 
terms are neglected in the k and e equations. Along the 
stagnation-point streamline, the following forms of the k and 
e equations are found: 


dk 2 k 2 d ( k 2 dk 

ay 7y-- 4aC »l + ‘-7y\ C ^dy, 


de 2 , e 2 d ( k 2 de\ 

ay—=- 4 a-c x cji + c 2 — - — c — — (6) 

dy k dy\ c 4 e dy) 

Introducing the coordinate transformations (eq. (1)) and the 
following forms of solutions for k and e, 

k = vaK(r) ) and e = va 2 E(r\) (7) 

results in the following k and e equations: 


dK K 2 d ( K 2 dK 

7 — = — 4 c u 1 - E c 

dr) E dij\ c 3 E dr j 


dE E 2 d ( K 2 dE\ 

+ (9) 


Equations (8) and (9) are solved for the turbulence properties 
K and E in the interval ijj > rj > rj 2 along the stagnation 
streamline, where rj ( is the boundary between flow regions 
1 and 2 and ij 2 is the boundary between flow regions 2 and 3. 


kr _ €f 

K = K, = — and E = E X = — - 
va va 


where k F and tp are previously determined by equation (2). 

(2) By neglecting the diffusion term in equation (9) at -q u 
this equation gives the following additional boundary condition 
at 7 = 7F 


d E 1 ( £f . „ 

X = “ ^--40,0^, 


(3) It is also further assumed (ref. 13) that 


dK n 

— = 0 at 77 = rj 2 
dr) 


The finite difference method is used to solve the boundary 
value problem in equations (8) to (12). Appendix C con- 
tains a description of the numerical method. The relation 
7i = V F /(va) 05 is found from equations (1) and (4), and t ? 2 
is resolved in the process of the numerical calculation. 

The solutions of the turbulence properties in flow region 2, 
obtained in the manner described previously, provide the 
turbulence quantities K and E at 7 = 7 2 - These turbulence 
quantities are required for the analysis of the boundary layer 
flow region. 

Region 3 

A steady, two-dimensional, incompressible turbulent bound- 
ary layer flow is assumed in this near-wall region. The mean 
flow properties are described by the boundary layer continuity, 
momentum, and enthalpy equations. The momentum and 
enthalpy equations contain the turbulence double correlations 
mv and vt. These double correlations are related to the higher 
order turbulence correlations and the mean flow properties 
through the transport of the turbulence. A method to model 
the turbulence double correlations was developed by Donaldson 
et al. (refs. 16 and 17). This method is used to formulate a 
theoretical analysis of the turbulence in this flow region. The 
Reynolds stress transport equations are derived by following 
the theory in reference 16. Assumptions are made for the 
turbulence closure relations and a modeling length scale A. 
A summary of the theoretical formulation and assumptions is 
given in appendix B. 

For the present boundary layer flow analysis, the mean 
velocity, U e = ax, is assumed along the boundary layer edge 
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between regions 2 and 3. The mean flow properties and Turbulence equations: 
turbulence double correlations are also nondimensionalized as 


U T 

F = — , H = - 


U, 


R x = — 2 , 

U: 


_ v 

Ry ~wr 


R z = - 2 

U 2 


dR x dR x 4 RfdU' U e d ( n , dR x \ 

di dr, U e d£ v dr, \ dr, / 


+ ^ + 2S(20 0 ' 5 - 

dr, 2 dr, 


0.5 


+ 2| 


(F - 3 R x ) vR 

3 At/, 


s ur 


e = - 


vt 

UJe 


+ 2 SR X [ -2 


8F 2v dF dr, 2v 2 


di t/, 5 t, dx t/ e 2 A 2 


Introducing the previous dimensionless variables and the 
coordinate transformations (eq. (1)) into the mean flow and 
turbulence transport equations (see appendix B) results in the 
following forms of the flow-field governing equations: 

Continuity equation: 


where 


dF dV, 

2£ — + —t + F = 0 

di dr, 


(13) 


di ; dr, v dr, y dr, ) dr, 2 


+ 2 $ 


vR” 5 ( R - 3 R y ) 
A U~ e 3 


+ 2 £R y 


(16) 


. .dF 2v dF 2v dF i )r, 2v 2 

X 2— + + - 


di U e dr, U e dr, dx U 2 A 2 


* 7 7 \ -a/ r j 


U e dx U e 


Momentum equation: 


'>tr dF +v dF 2 «1-F 2 )dU e 

2 ^TZ + V t— = 


di dr, 


Ue a« 

d 2 F „ n , as 
+ -^ + (20 0 - 5 - 
or? 07? 


(14) 


(17) 


„ 3F, „ dR z 2 R z FdU e 

2 ZF—1 + V . — - = - 2f — — — + 

di dr, u e di 


¥l±fx R 0-5 d A\ 

v ai/V dr,) 


+ 


d z R, 


vR 05 (R - 3 R z ) 


+ 2$ 


dr, 2 ’ AU e 
2v 2 

-2iR z ^~. 

Z U 2 A 2 


(18) 


Enthalpy equation: 


„dH dH 

2 %F — +V t — = 
d^ dr, 


2 jFH dT e 
T e di 


+ OS* 

Pr dr, 2 


+ (2 i) 05 — 

orj 


(15) 


Slfaui?! 

du dr, U e di v dr, \ dr,/ 


d 2 S n , dF 

+ 5 ? + CS) *’ 7 , 


+ 2(Sl - 


vR 05 2v 2 


UtA U 2 A 2 


(19) 
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and 


Similarly, 


Fd dT e U e d in0S dd 

- 2£ + 3 — — A R 05 — 

T e d£ v dr, \ dr, 


dR, 

~ = R zf\ ( Vo) 

dk 


U U n c (7/7 

+ -^ + (20°^.— 

07 ] 07 ] 

l 2 i -e( aF i v dFdri vR ° 5 21,2 
u e dr, 8x U e A C/ 2 A 2 


dS 

— — Sf 2 ( <p 0 ) 
dk 


where R = R x + R y + R z . 

In order to solve the previous governing equations for the 
mean flow and turbulence properties, the initial and boundary 
conditions are derived analytically in the following manner. 

Initial conditions.— The initial profiles of the mean flow 
properties and the turbulence correlations at <p = <po are 
determined analytically from equations (13) to (20). The 
following assumptions are made: 

(1) At 1,2 location, 


aF = o 

dr, ’ dr, 

dR x dR dR z 

— - = 0, — = 0, — = 0 

dr, dr, dr, 

dS „ dd 

— = 0, — = 0 

dr, dr, 


dd 

— = 0M<fio) 
dk 


are also assumed at <p 0 . 

(3) The mean flow is symmetric about the <p = 0 plane. In 
order to facilitate the analysis, this mean flow behavior is also 
assumed at small distances away from the stagnation point, 
<p = 0. Therefore, dF/d£ = 0 and d///d£ = 0 are imposed 
along the r, direction at the <p 0 location. 

(4) The conditions F = 0, H = TJT e , R x = 0, R y = 0, 
R z = 0, S — 0, and d = 0 are also imposed at the cylinder 
surface. 

By applying these assumptions (eqs. (21) to (24)) to the 
turbulence governing equations, the following relations are 
found at <p = <p 0 and rj = rj 2 : 


^x.o.e 


3 U 0 , e \ e 


-2 | 4(dU/dQ 0 


H — q 

u 0 Ae ulX 


In addition, 


F= 1, H= 1, S = c 5 (R x R y ) 05 , and d = 0 

are imposed at r, 2 . Also, c 5 is an empirical constant; the value 

c 5 = 0.001 is used in this study. 

(2) The Reynolds normal stress components u 2 and v 2 are 
assumed to be independent of ip for very small <p. Therefore, 
at <p 0 . 


D 

t'y.Oye 


3U oe A e 


Re<Po V oe A e U\ e A 2 


, , , 2 (dUm o .e , 

/i(^o) = + 


3^ 0 .A, 0 ,A 


flip o) — 


2 (ftV/dS )" _ jKl _ 2» 2 

Uo, c u 0 ' A ulAl 


y \Revl 


where R o e = 2mK 2 /Ul e is assumed. Since d = 0 is imposed 
at r, 2 , dd/d£ = 0 occurs at the boundary layer edge. This zero 
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(36) 


gradient condition is also used within the boundary layer at 
<p = (p 0 . It follows from the assumption (eq. (25)) that 

Mv o) = 0 (30) 


With the aid offi(<p 0 ),f 2 (<p 0 ), and / 3 (^ 0 ). assumption 2 defines 
the mean flow, streamwise turbulence gradient terms. 
Substituting the mean flow, streamwise gradients into 
equations (13) to (20) results in a set of ordinary differential 
equations with the boundary conditions specified at i; = 0 and 
i] = i/ 2 . These equations are solved numerically for the 
boundary layer flow property profiles at y> 0 = 0.04° in order 
to avoid the singularity at the stagnation point {<p = 0°). These 
property profiles are used as the initial conditions for the down- 
stream, boundary layer flow-field analysis. 

Boundary conditions —The turbulence double correlations 
along the boundary layer edge are also required to analyze 
the flow field. The turbulence kinetic energy is first defined 
by using the k-e two-equation turbulence modeling. The 
turbulence kinetic energy is then substituted into the Reynolds 
stress equations to determine the turbulence double 
correlations. 

By following the k-e two-equation turbulence modeling 
technique (appendix A), the turbulence modeling equations 
can be written in the x-y coordinate system as 


dk 2 k 2 
ax — = 4 a c.. — 
dx * e 


d 

e + — 
dx 


( k 2 dk\ 
Vc 3 edx) 


(31) 


de 
ax — 
dx 


= 4 a 2 C\C fl k 


e 2 d ( k 2 de\ 
C2 T + ^\ ' l ~^edx) 


(32) 


If the coordinate transformations (eq. (1)) and the following 
forms of solutions for k and e 

k = vaK(i) and e = va 2 E(£) (33) 


are introduced into equations (31) and (32), the turbulence 
modeling equations become 


^dK K 2 „ d 

2£ — = 4c E + 2£ — 

dk * E dk 


\ c 3 E dk) 


(34) 


dE E 2 d 

V -,4c,c,K-c 1J + V - 


C *c 4 Edk) 


(35) 


The analysis of the turbulence in flow region 2 defines the 
following initial conditions: 


K = K 2 and E = E 2 

at k = £2 an d V = r) 2 . 

With the assumptions that the gradients dRJdi j, dRyldr), 
and dR z /dr\ vanish at the boundary layer edge, equations (16), 
(17), and (18) give the following theoretical relations among 
the mean velocity, the Reynolds normal stresses, and their 
mean flow streamwise gradients: 


dR 


1.5 


x,e _ 


dk 


4R x , e dU e | v_ R l 
U e dk u e 3A e 


vR^Rr, 2v 2 


U e K 


u 2 Ar x ’ e 


(37) 


dR, 


dk 


y>e _ 


vR' e 5 v R°%, e 2v 2 „ 

: A, 


3U e A e U e A e 


2 * 2 y, e 


Ui A e 


(38) 


and 


dR. 


•Z.e _ 


dk 


2R z e dU e + vR, 


1.5 


U e dk 31/A 


pR 


0.5 


UAe 


-R, , — 


2v L 


;R- 


u 2 a 2 z ’ e 


(39) 


The previous equations can be solved for the Reynolds 
normal stresses. The sum of the Reynolds normal stresses R 
should be consistent with the turbulence kinetic energy 
obtained from equations (34) and (35). The Reynolds normal 
stresses along the boundary layer edge are determined from 
this consistency requirement. Therefore, the following analysis 
is used. 

The previous analysis of the initial profiles defines the 
Reynolds normal stresses at k - £2 an d V = rj 2 . These are the 
initial conditions to integrate equations (37) to (39) along the 
downstream direction. In this study, the fourth-order Runge- 
Kutta numerical method is used. The turbulence kinetic energy 
obtained previously from equations (34) and (35) is used to 
evaluate the derivatives. The integration step is adjusted until 
the turbulence kinetic energy (R xe + R ye + R ze )U 2 /2 
satisfactorily converges to the corresponding value predicted 
by the k-e two-equation turbulence modeling. The cor- 
responding Reynolds normal stresses are the boundary 
conditions at the boundary layer edge. 

Finally, the following conditions are also imposed: 

(1) At the boundary layer edge 77 = r/ 2 , F e — 1, H e = 1, 
S e = c 5 (R xe R y( ) 05 , and 6 = 0. 

(2) At the surface tj = 0, F = 0, H = TJT e , R x = 0, 
R y = 0, R. = 0, S = 0, and 6 = 0. 
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A numerical computational procedure to calculate the 
boundary layer flow field was developed in a previous study 
(ref. 18). This numerical method and the other computational 
procedures used in this study are described in appendix C. 


Results and Comparison with Experiment 

In order to verify the previous theoretical analysis, a 
numerical computation was performed to predict the existing 
measurements. A set of detailed experimental data on the effect 
of free-stream turbulence on heat transfer from heated 
cylinders placed normal to the airstream was reported by 
Lowery and Vachon (ref. 5). Tests were conducted over a 
range of Reynolds numbers and at different turbulence levels. 
The two-point turbulence longitudinal correlations in the free 
stream were measured. These data can be used to determine 
the turbulence longitudinal microlength scale that is required 
to start the present numerical calculations. Thus, the test 
conditions of this reference experiment were used as input to 
the present numerical computational procedures. 

Examples of the free-stream turbulence longitudinal 
correlation and the corresponding turbulence longitudinal 
microlength scales, with Tu = 0.013, 0.017, and 0.028, and 
Re = 2.50x10 s , are shown in figure 2. The microlength 
scales were obtained by curve-fitting the experimental data 
with equation (3). On the basis of the free-stream mean veloc- 
ity, the free-stream turbulence intensity, and the turbulence 
longitudinal microlength scale, the present numerical 
computational procedures were performed to calculate the 
turbulence kinetic energy, the turbulence double correlations, 



Figure 2.— Free-stream turbulence longitudinal microlength scale. Reynolds 
number Re, 2.50 x 10 5 . 


and the mean flow properties within the stagnation flow field. 
The surface heat transfer rate was also calculated from the 
mean temperature distribution. In the following sections, these 
numerical results are described and are compared with the 
experiment of reference 5. 

Turbulence 

To show the variation of the turbulence kinetic energy along 
the stagnation-point streamline, some analytical results are 
plotted in figure 3. Except for the lowest turbulence case, the 
turbulence kinetic energy increased within the inviscid flow 
region in the direction toward the surface and reached a 
maximum value at the boundary layer edge. The turbulence 
kinetic energy then decreased continuously and vanished at 
the cylinder surface. This turbulence kinetic energy variation 
is similar to the result in reference 13. However, the present 
analysis predicts a reduction in the turbulence kinetic energy 
along the stagnation-point streamline when the free-stream 
turbulence intensity and Reynolds number are small 
(Tu = 0.013 and Re = 1.90x10 s ). 

Examples of the Reynolds normal stress profiles within the 
boundary layer at the stagnation point are shown in figure 4. 
These results are for the cases with Re = 2.50x 10 5 and 
Tu = 0.017 and 0.028. The Reynolds normal stresses 
decreased continuously from the boundary layer edge in a 
direction normal to the surface. Large values of R y among the 
Reynolds normal stresses were calculated within a small region 
near the surface. These Reynolds normal stress profiles 
differed significantly from those of a flat-plate turbulent 
boundary layer (ref. 20), where the largest Reynolds normal 
stress was measured in the x direction. 



Figure 3.— Turbulence kinetic energy along stagnation streamline. 
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MEAN FLOW STREAMWISE VELOCITY, F = U/U e 


Figure 4.— Reynolds normal stress profiles within stagnation-point boundary 
layer. Reynolds number Re, 2.50x10 s . 


Figure 5.— Mean flow streamwise velocity component within stagnation-point 
boundary layer. Reynolds number Re, 2.50 x10 s . 


The analytical predictions of the mean velocity component, 
F = UIU e , within the boundary layer at <p = y> 0 are shown in 
figure 5. These profiles are compared with the mean velocity 
profile of the theoretical laminar boundary layer. The 
comparison shows that the free-stream turbulence increased 
the velocity near the surface and decreased the velocity near 
the boundary layer edge. Existing measurements (ref. 21) of 
this velocity component with Re = 2.50x 10 5 and Tu ~ 0.05 
are also plotted in this figure. The analytical results are in good 
agreement with the measurements. 

Examples of the computational results of the turbulence 
kinetic energy and the Reynolds normal stresses along the 
boundary layer edge are shown in figure 6. In this figure, the 
turbulence properties are nondimensionalized with the corre- 
sponding stagnation-point values. The turbulence kinetic 
energy increased continuously along the downstream direction. 
The rate of change in the kinetic energy is a function of the 
free-stream turbulence properties \ F , Tu, and Re. A large 
increase in the turbulence kinetic energy was found for the 
case with small free-stream turbulence intensity and a small 
Reynolds number. The Reynolds normal stress components 
Ug and Vg increased at the downstream locations. However, 
the Wg component was reduced near the stagnation point. The 
effect of the microlength scale on the Reynolds normal stress 
along the boundary layer edge is shown in figure 7. The 
computational results for the cases with Re = 2.50 x 10 5 and 
Tu = 0.028 but different length scales, \f/D = 0.0053 and 
0.0079, are presented. The turbulence microlength scale 
affected the Reynolds stress dist ribut ions. Increasing the length 
scale decreased the ratio of Ug/ul e and increased the ratio of 


Wglwlg along the boundary layer edge. No significant 
difference in the ratio of vj/vl e was calculated. 

Surface Heat Transfer 

The present analytical results of the surface heat transfer 
are defined in terms of the Nusselt number 



(40) 


The mean temperature gradient ( dTldy) y=0 at the surface was 
calculated numerically from the mean temperature profile, 
which in turn, was obtained previously from the solution of 
the boundary layer flow. 

The analytical results for the stagnation-point surface heat 
transfer rate over a range of the free-stream Reynolds number 
and the turbulence intensity are compared with the 
experimental data (ref. 5) and existing correlations (refs. 5 
and 6) in figure 8. The surface heat transfer rate is expressed 
in the form of the Froessling number (Nu Re~ 0 5 ). Linear 
variation in the turbulence longitudinal microlength scale as 
a function of the free-stream Reynolds number was also 
assumed in the present calculation. The turbulence longitudinal 
microlength scale decreased as the Reynolds number increased. 
With constant free-stream turbulence intensity, turbulence 
microlength scales within ± 10 percent of the values in figure 2 
were used in the numerical computation for a range of the free- 


9 


TURBULENCE KINETIC ENERGY. K./K 


REYNOLDS 

DIMENSIONLESS 

DIMENSIONLESS 

DIMENSIONLESS 

INITIAL BOUNDARY LAYER EDGE 

NUMBER, 

TURBULENCE 

FREE-STREAM 

FREE-STREAM 

STATION 

Re 

LENGTH SCALE, 

TURBULENCE 

TURBULENCE 




Ap/D 

KINETIC 

DISSIPATION 

DIMENSIONLESS 

DIMENSIONLESS 



ENERGY, 

RATE, 

TURBULENCE 

TURBULENCE 



k f 

e f 

KINETIC 

DISSIPATION 





ENERGY, 

RATE, 





K o,e 

E o,e 

1.96X10 5 

0.0053 

12.27 

11.45 

7.51 

3.79 

3.00 

.0053 

19.06 

11.45 

19.03 

11.42 

/ 




/ 


/ 

^K*/K 0 .e 


V ul/ujj.e 






>-ve^vo.e 

r K e /K o,e 

l 


1.96XKP 0.0086 

2.50 .0053 


sf~ u e /u o.e 


.04 . 06 . 08 .1 


.4 .6 ,8 1 2 

AZIMUTH ANGLE. <P . DEG 


(a) Turbulence intensity Tu, 0.013. 

(b) Turbulence intensity Tu, 0.028. 

Figure 6. — Turbulence kinetic energy and Reynolds normal stresses along boundary layer edge. 




Figure 7.— Effect of turbulence longitudinal microlength scale on Reynolds normal stresses along boundary layer edge. Reynolds number Re, 2.50 x10 s ; 
free-stream turbulence intensity Tu, 0.028. 



Figure 8.— Comparison of stagnation-point Froessling number between 
analysis and experimental correlations. 


stream Reynolds number of 1.50 </?eXl0 -5 < 3.00. 
Because it includes the effect of the free-stream turbulence 
microlength scale, the present analysis calculates a larger 
stagnation-point surface heat transfer rate than is indicated by 
the correlations in references 5 and 6. 

Examples of the computational results of the turbulence 
double correlation —vtlU e T e within the boundary layer at the 
stagnation point are shown in figure 9. The maximum 
turbulence correlation occurs near the surface. Increasing the 
free-stream turbulence intensity also increases the turbulence 
correlation vt. Together with the characteristics of the Reynolds 
normal stress component (fig. 4), the present analysis indicates 
that the free-stream turbulence penetrates into the near-surface 
flow region and induces a large turbulence correlation vt. Thus, 
a surface heat transfer rate larger than that for laminar flow 
occurs. 

Some analytical results of the surface heat transfer rate 
within a small distance from the stagnation point are plotted 
in figure 10 in terms of the Froessling number. An 
approximately constant Froessling number is calculated only 
for the case with Tu = 0.013. For Tu = 0.017 and 0.028, the 
analytical results show a reduction in the Froessling number 
near the stagnation point. With high free-stream turbulence 
intensity, a large turbulence dissipation rate may occur near 
the stagnation point. This turbulence dissipation rate reduces 
the rate of the increase in the turbulence kinetic energy. This 
rate of increase in K e is not enough to induce a constant 
surface heat transfer rate. 

To investigate the effect of the turbulence longitudinal 
microlength scale on the surface heat transfer rate, we 


li 



Figure 9.— Profiles of turbulence correlation vt within stagnation-point 
boundary layer. 


compared the analytical results of the turbulence kinetic energy 
and the Froessling number for \ F ID = 0.0053 and 0.0079 
with Re = 2.50x 10 5 and 7m = 0.028 in figure 11. With 
\ f ID = 0.0079, the present analysis computed the 
experimental data in reference 5. A lower level of the 
turbulence kinetic energy and a smaller Froessling number 
were also predicted for the case with a small free-stream 
turbulence longitudinal microlength scale. These computational 
results also show a reduction in the Froessling number near 
the stagnation point. 

In this study, the approximation 

S e = c 5 (R Xie R y , e ) 0 - 5 

was used along the boundary layer edge. This approximation 
is similar to the correlation observed from the flat-plate 
turbulent boundary layer study (ref. 20). In the present 
analysis, c 5 = 0.001 was found to give the best computational 
results. In an earlier study (ref. 18), a different value of c 5 
(0.0001) was used because the boundaiy layer edge location 
and the corresponding turbulence kinetic energy (or the 
Reynolds normal stresses) were not analyzed theoretically. If 
one uses the smaller value of c 5 , then a small (nonzero) nv 
is imposed along the boundary layer edge. This approximation 
is different from the k-e two-equation turbulence modeling, 
which requires S e = 0 for the external flow field in this study. 
However, the requirement of S e = 0 is itself an assumption. 

The assumption vt e = 0 was also imposed at the boundary 
layer edge. Intuitively, the turbulence dynamic equations can 
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Figure 10.— Boundary layer edge turbulence kinetic energy and corresponding Froessling number. 
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Figure 11. — Effect of turbulence longitudinal microlength scale on turbulence kinetic energy at boundary layer edge and Froessling number. Reynolds 
number Re, 2.50x10 s ; turbulence intensity Tu, 0.028. 


also be used to study the variation of the turbulence double 
correlation vf within the external flow region. If vf is nonzero 
in the external flow field, then a different vf distribution may 
occur within the boundary layer flow region, and the present 
analysis will predict different stagnation-flow surface heat 
transfer rates. If this is true, it might resolve the different levels 
of the stagnation-flow surface heat transfer rate reported in 
various experiments. Although the task requires the vt 
boundary condition within the free stream, the effect of the 
free-stream temperature fluctuation on the surface heat transfer 
rate of the stagnation flow can be studied by using the present 
analysis. 

Conclusions 

A theoretical analysis was formulated to model the effect 
of the free-stream turbulence on the surface heat transfer of 
a stagnation-flow region. The flow field of interest was around 
the forward stagnation point of a circular cylinder in a turbulent 
flow. The k- e two-equation turbulence modeling was used to 
analyze the turbulence kinetic energy within the external 
inviscid flow region. The time-averaged turbulence transport 
equations were employed to predict the turbulence double 
correlations within the boundary layer. From the boundary 
layer edge conditions of the turbulence kinetic energy, the 
transport equations of the mean flow and turbulence were 
solved numerically for the turbulence double correlations, the 
mean velocity, and the mean temperature. 


The present analytical results show that, depending on the 
free-stream turbulence properties, different levels of turbulence 
kinetic energy occur along the boundary layer edge. These 
different levels of turbulence kinetic energy induce different 
surface heat transfer rates. Increasing the free-stream 
turbulence intensity or the turbulence longitudinal microlength 
scale also increases the stagnation-region surface heat transfer 
rate. It is postulated that the free-stream turbulence penetrates 
into the boundary layer and induces high surface heat transfer 
rates. 

The analysis presents a procedure to calculate the stagnation- 
point surface heat transfer rate as a function of the free-stream 
turbulence intensity, the turbulence longitudinal microlength 
scale, and the Reynolds number. Because of limited data on 
the turbulence longitudinal microlength scale, this analysis was 
verified only for some cases with a Reynolds number of the 
order of 10 5 and a free-stream turbulence intensity of less 
than 5 percent. By considering the free-stream turbulence 
longitudinal microlength scale, this analysis predicts a higher 
stagnation-point surface heat transfer rate than was obtained 
from the existing correlations relating the free-stream 
turbulence intensity, the Reynolds number, and the surface 
heat transfer rate. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, June 2, 1987 
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Appendix A 

k-e Two-Equation Turbulence Modeling 


Launder and Spalding (ref. 22) developed the k-e two- 
equation turbulence modeling technique. It has been found 
effective by many researchers in modeling the turbulence of 
different turbulent flow configurations. In the present study, 
this technique is also used to model the turbulence along the 
stagnation-point streamline and the boundary layer edge. 


k-e Transport Equations 

With a two-dimensional mean flow field, the modeling 
equations for the turbulence kinetic energy k and the turbulence 
dissipation rate e can be written in an x-y orthogonal coordinate 
system as 


jjdk + v <tt = _a/^dA + _a_A^'\ 

dx dy dx \c 3 dx) dy \c 3 by) 


-^dU 
— u — 
dx 


analysis: c 3 - 1.44, c 2 = 1.92, c 3 = 1.00, c 4 = 1.30, and 
= 0.09. 

Applications of k-e Equations 

The k and e equations are used to predict the turbulence 
kinetic energy along the stagnation-point streamline in flow 
region 2 and the boundary layer edge between flow regions 2 
and 3. The mathematical formulations of the k and e equations 
are described in the following section. 

Along Stagnation-Point Streamline 

The turbulence kinetic energy variation along the stagnation- 
point streamline is analyzed. For the stagnation-point flow 
field, the inviscid flow analysis predicts the mean flow velocity 
components as 


- MV 


/ dV dU\ 
dy) 



(Al) 


U = ax and V = — ay (A4) 

with a = 4V f /D. Therefore, the velocity gradients 


and 


u* + v?i 

dx dy 


d / v, de\ + d / v, d«\ 
dx \c 4 dx) dy \c 4 dy) 


e^dU e — 

— Ci - M Ci - MV 

k dx k 


fdV dU\ 

\3jc dy) 



dU n dV 

— = 0 , — 
dy dx 



(A5) 


are found from equation (A4). 

Substituting the previous relations into equations (Al) and 
(A2) results in these equations being simplified and becoming 


e-^dV e 2 

■Cl - v c 2 — 

k dy k 


(A2) 


dk , k 2 

ay— = -4 a c — + e 
dy e 


d k 2 dk 

— 1 c i* 

dy\ c 3 e dy 


(A6) 


The turbulent dynamic viscosity v, is related to k and e by 
v, = c^ 2 /e. The closure assumptions for the Reynolds 
stresses are 


de 


ay— = - 4 a 2 c l c )l k + c 2 — - — — 


dy 


k 2 de 


k dy\ c 4 e dy 


(A7) 




-5- „ dU 

2 


— k 

dx 

3 

— (dU 


MV = V,\ 

+ — 

'\dy 

dxj 


-t , dV 2, 

- v 2 = 2v, k 

dy 3 


J 


(A3) 


and c v c 2 , c 3 , c 4 , and c M are empirical constants. The 
following values are usually used in the turbulence flow 


Furthermore, by introducing the coordinate transformations 
(eq. (1)) into the previous equations, the k-e two-equation 
turbulence modeling becomes 


arj 


dk 

dr) 


= — 4 a 2 c.. 


a d 

+ e 

v dr) 


( k 2 dk\ 


(A8) 


de . , , e 2 

ar)—=~ 4a C\C k + c 2 — 
dr) k 


ad / k 2 de 
v dr) \ 11 c 4 e drj 


(A9) 


14 



We are seeking the solutions of k and e in the forms 

k = vaK{-q) and e = va 2 E{ rj) (A10) 

Substituting the previous definitions into equations (A8) and 
(A9) results in the following dimensionless turbulence 
modeling equations being obtained along the stagnation-point 
streamline: 


The closure assumptions (eqs. (A3)) are also used in deriving 
the previous k and e equations. The mean velocity properties 
in equations (A 13) and (A 14) are substituted with the results 
in equations (A4) and (A5). The k and e equations become 



dK 


V 


dr ) 



dK\ 
c 3 E dr] J 


(All) 


dE E 2 d ( K 2 dE 

V — = ~ 4ciC„a + c 2 c„ 

drj * K dr]\*c 4 Edr] 


(A 12) 


Along Boundary Layer Edge 

Equations (Al) and (A2) are also used to analyze the 
turbulence kinetic energy along the boundary layer edge. The 
basic assumptions are that the gradients dkldy and del by vanish 
at the boundary layer edge. Thus, the k and e equations become 


T dk k 2 ( dU \ 2 fdvV fdU dV \ 2 

U— = c — 21 — +2 — + — + — 

dx * e \dx J \dyj \dy dx J 

d_( k^dk\ 
dx \ M c 3 e dx ) 


(A13) 



e 2 d ( k 2 3e\ 
C2 ~k + !bc 


(A 14) 


de 


ax — = 4a 2 c,c.,k — c 7 1 ( c„ 

dx M k dx\*c 4 edx 


k 2 de 


(A 16) 


Again, the x coordinate is transformed by using equation (1). 
The following forms of solutions for k and e are also sought: 


*(£) = vaK{%) and e(|) = va 2 E(^) 


(A17) 


Paralleling the results in flow region 2, the following 
dimensionless forms of the k and e equations are obtained 
through mathematical manipulations: 


. y dK M K 2 „ „ d( K 2 dK 

2£ — = 4c E + 2£— ( c„ 

dH » E dA*c 3 Ed% 


(A 18) 


. y dE E 2 d( K 2 dE\ 

2 ( -, ic ,c f K-c 2J + 2 ( -U--\ «A,9, 
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Appendix B 

Turbulence Transport Equations 


For engineering applications, the concepts of eddy viscosity 
and eddy conductivity have been used by researchers to 
calculate the momentum and heat transfer in turbulent flow. 
However, in cases where the dynamics of turbulence is 
important, the concept of eddy transport may not be sufficient 
to study the turbulence. Thus, Donaldson et al. (refs. 16 and 
17) developed a method, the method of invariant modeling, 
to model the dynamics of the turbulence correlations. It is used 
in this study to analyze the turbulence within the boundary 
layer flow region. 


and 


d — d _ 

U— vt + V— vt - 
dx dy 


-jdT _ dU 

v V vt — 

dy dx 


1 d — 1 dt d 

pt + - p vvf 

p dy P dy dy 


d 2 — dv dt 

+ v — r vt - 2v 

dy 2 dx, dx, 


(B5) 


Turbulence Transport Equations 

The one-point turbulence dynamic equations are used in the 
method of invariant modeling. The details in reference 16 are 
followed to derive the transport equations for the turbulence 
double correlations u 2 , v 2 , w 2 , uv, and vt in the x-y-z 
coordinate system. These equations are 


dx dy 


— dU — j dU 2 du 
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dy dx p dx 


d d 2 u 2 du du 

VMM + V r 2v 
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(Bl) 


dv 2 
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U V V — = 2 v 2 1 — p vvv 
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2d — d 2 v 2 dv dv 

-I VD + V ; 2v 
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(B2) 


Turbulence Closure Assumptions 

Assumptions for the higher-order turbulence correlation 
terms in equations (Bl) to (B5) must be made so that these 
turbulence equations and the mean flow transport equations 
can form a closed system of equations to describe the flow 
field. These closure assumptions are based on the turbulence 
double correlations and a length scale A. In this study, the 
following forms of the closure assumptions (refs. 16 and 17) 
are used: 
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Figure 12.— Modeling length scale of turbulence closure assumption. 


constant. This length scale profile is shown schematically in 
figure 12. The pivot point location y p is equal to 10 percent 
of the theoretical laminar stagnation-point boundary layer 
thickness 6. At the boundary layer edge, we found that the 
length scale A f can be related to the free-stream turbulence 
intensity by 

A^ = (0.25 - 3Tu)8 (B6) 


bu du u 2 bv bv v 2 

bXi bx, A 2 bxj bxj A 2 


bw bw w 2 

bxj bx, A 2 ’ 


bu bv uv 

bXi bxj A 2 


Applications of Turbulence 
Transport Equations 

In the present study, the mean flow is two dimensional, 
steady, and incompressible. It is described by the following 
turbulent boundary layer continuity, momentum, and enthalpy 
equations: 


dv bt vt 

bXi bx j A 2 


The length scale distribution within the boundary layer in 
flow region 3 of this study is also assumed to be similar to 
the mixing length profile of a flat-plate turbulent boundary 
layer. Within a small distance from the surface, the length scale 
varies linearly as a function of the normal distance from the 
surface. Away from the surface, the length scale remains 
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The turbulence double correlations uv and vt are modeled 
by using the turbulence dynamic equations. They are obtained 
by substituting the turbulence closure assumptions into 
equations (Bl) to (B5). Thereby, the turbulence dynamic 
equations are found to be 
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Together with the length scale assumptions, the previous 
equations form a closed system of turbulence transport 
equations. These equations are used to analyze the turbulence 
and the mean flow properties in the boundary layer flow 
region. They are nondimensionalized according to the 
procedures in reference 23. The resulting dimensionless 
equations were reported in a previous study (ref. 18). 
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Appendix C 

Numerical Computation 


The finite difference method is used to calculate the 
turbulence and the mean flow properties within the flow field. 
Summaries of the finite difference scheme and the numerical 
computational procedures are presented in this appendix. 

Finite Difference Scheme 


The following central finite difference scheme with variable 
grid size in the direction is used. The derivatives at each 
grid point m are approximated by 


G m + 1 — 


Vm + 1 Vm — 1 


All of the transport equations (eqs. (All), (A 12), (A 18), 
(A 19), and (14) to (20)) can be represented by the following 
equation: 


dG 3G d ( dG 

A G 1 ~Z + A G1 A CiG + A g 4 — 1 A G5 — 
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where G is either a turbulence or a mean flow property. A 
comparison between equation (Cl) and each transport equation 
shows that A a (i = 1 to 9) may be constant, or the functions 
of the properties and their gradients, or the coordinates £ and 
r). The transport equations are written in the forms of equation 
(Cl) to describe the finite difference scheme. 

k-e Turbulence Modeling Equations 

In flow region 2, either of the dimensionless equations (All) 
and (A 12) can be written as 
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with G = K or E. 

Similarly, equations (A 18) and (A 19) in flow region 3 can 
also be written as 


Substituting the previous expressions into equation (C2) 
yields a finite difference equation. 


dG „ d dG\ 

A G1 — - A G3 G + A G4 ~ I Aq 5 — 1 + Acg 

d£ d£ \ dU 


with G = K or E. 

Equation (C3) has a form similar to equation (C2) regardless 
of the difference between the £ and ?/ coordinates. A similar 
finite difference scheme is used to convert equations (C2) and 
(C3) into finite difference forms. As an illustration, the steps 
in transforming equation (C2) are described in the following 
paragraphs. 


BiGm+i + G m + 5 3 G m _ i - B 4 (C6) 

at each grid point m. The coefficients S, (i — 1 to 4) are 
functions of the various coefficients A a and the grid size at 
the grid points m + 1, m, and m — 1. 

Applying the previous mathematical manipulation at every 
grid point within the interval rji > rj > r) 2 results in a set of 
finite difference equations for equation (C2). The boundary 
conditions of G are specified at i)\ and rj 2 . This set of 
equations is solved for the G values at each grid point m. 
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Boundary Layer Flow Transport Equations 

All of the governing equations (eqs. (14) to (20)) can be 
written in the following form: 


Substituting the previous relations into equation (C7) 
produces a set of finite difference equations, similar to the form 
of equation (C6), for each of the governing equations (eqs. 
(14) to (20)). 


dG dG d 

A G1 — + A G2 — + ^G6 ~ 

a £ or) or) 



Computational Procedure 


d^G 

+ A CS ~r~2 + (C7) 

01) 

where G can be one of the mean flow properties or the 
turbulence double correlations. 

For the initial profile calculation, the assumptions and the 
results in the initial profile analysis are used to define 3G/d£. 
For the downstream calculations, the gradient 3G/d| is 
approximated by 


(3G\ = ( % -Gg) 

A? 

where G £ is the G m at the previous station. Constant grid size 
is used in the i) direction. The derivatives in the ij direction 
are approximated with the following central finite difference 
scheme: 


dG \ _ G m+ i - <j m _i 
di)) m 2A i) 


(C9) 


( d 2 G \ _ G m+ 1 - 2 G m + G m _! 

\ d V 2 /m ( A, j) 2 


and 


dG 

d ( 4,77 — 
01 ) 


di) 


(CIO) 


[4g7., 


+ \/2'~ r m+\ ' 


(4 g 7 


m + 1/2 


+ A ( 


G1 ,m— 1/2 


) G m A q im— 1 llG m _ i ] 


(Ar ,) 2 


(C11) 


The sets of finite difference equations are solved by using 
iterative numerical procedures. Depending on the boundary 
conditions, slightly different iterative schemes are employed 
in the numerical computations. 

k-e Turbulence Modeling Equations 

Different procedures are used in flow regions 2 and 3. They 
are described separately in the following sections. 

Flow region 2.— Equations (10) to (12) describe the 
boundary condition at i) = i)\ and i) 2 . However, i) 2 is not 
defined. Therefore, t) 2 is first assigned an arbitrary value to 
start the numerical calculation. The existing analysis (ref. 13) 
indicates that t) 2 is approximately 5. Linear variations in K 
and E are first assumed. The coefficients of the finite difference 
equations are evaluated with these K and E assumptions. Thus, 
two sets of linear algebraic equations are obtained to relate 
K and E at the grid point locations. Since K and E are 
prescribed at i) X and i) 2 , these equations are solved for K and 
E distributions by using the method of successive substitution. 

The previous procedures are repeated until K and E converge 
at each grid point. The boundary values of K and E at i) 2 are 
adjusted to ensure that the K and E profiles satisfy the boundary 
conditions (eqs. (11) and (12)). 

The convergence in K and E is sensitive to their assumed 
values at the start of the computations. On the basis of k-e 
two-equation turbulence modeling, Strahle (ref. 15) developed 
analytical links between the turbulence properties k and e at 
the stagnation point and in the free stream. By neglecting the 
diffusion terms and using the same value for the empirical 
constants c x and c 2 (i.e. , c x = c 2 = 1.90), these analytical 
links can be written as 



(C12) 


and 


with 


’G7,m+ 1/2 


'C7.B+ 1 


+ A, 


Gl,m 


e 2 - 


2.4 V F k F 
D 


2.111 


- 1.111 

e F 


(Cl 3) 


4(77, m— 1/2 


iGH.m 


+ A, 


GJ ,m — 1 


2 


According to equation (7), the previous relations are 
nondimensionalized and they become 
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K 2 = K} m (C14) 

and 

E 2 = (0.6K f ) 7 (C15) 

With the isotropic assumptions (eq. (2)) in flow region 1, 
the approximate boundary conditions for K and E at i) 2 are 
calculated from the previous relations. These K and E values 
are used to start the numerical computations. 

The grid size within the computational interval is generated 
by using the model in reference 23. Three hundred to four 
hundred grid points are used in the computations. The grid 
is tightened near tj, and rj 2 to ensure the boundary conditions. 

Flow region 3.— The aforementioned numerical calculation 
provides the conditions of K and E at £ 2 . From these 
conditions, the K and E along the boundary layer edge are 
also calculated by using the sets of finite difference equations 
corresponding to equations (34) and (35). The computational 
domain £ 2 < £ < £ 3 is limited to a small region near the 
stagnation point. In this study, £ 3 , corresponding to <p = 15°, 
was chosen. 

The numerical procedure is similar to that used in flow 
region 2. First, K and E at the £ 3 location are estimated and 
linear distributions in K and E between £ 2 and £ 3 are 
assumed. The coefficients of the finite difference equations 
are calculated by using the assumed K and E distributions. Two 
sets of linear algebraic equations, relating the K and E at the 
grid points, are obtained. These algebraic equations are solved 


with the method of successive substitution to update the K and 
E distributions. This iterative computational procedure is 
repeated until K and E converge at each grid point. The 
numerical values of K and E at £ 3 are also adjusted in each 
iteration to facilitate the convergence in the numerical 
computations. 

The grid size within the computational interval is also 
generated with the technique used in flow region 2. Two 
hundred to three hundred grid points are used in the present 
computations. The grid is also tightened near £ 2 and £ 3 . 

Boundary Layer Flow Equations 

The numerical computational procedure was developed 
previously in reference 18. It contains two main routines at 
each £ downstream location. First, the velocity field is 
calculated. This velocity field is then substituted into the energy 
equation to compute the temperature distributions. 

Similar routines are repeated at the next downstream 
locations. At each £ station, the property G at the previous 
£ station is used to define the coefficients A Ci in the first 
iteration, and the updated value for G at the current station 
is then used to determine these coefficients in the next iteration. 
However, in the initial profile calculations, the coefficients 
A a are calculated by using the laminar stagnation-point 
boundary layer flow profiles. 

For the cases calculated in this study, the computational 
domain is limited to 0.04° < tp < 15° and i? < 6.0. Constant 
grid size (Aij = 0.05) in the ij direction and variable step size 
along the £ direction (A£ corresponding to A <p = 0.01° for 
ip < 0.2° and A <p = 0.2° for <p > 0.2°) are used. 


\ 
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